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Abstract 

In  this  paper,  a  comparison  between  the  Entropy  Generation  Minimization  method  and  the  Power 
Maximization  technique  is  presented.  The  assessment  is  performed  by  analyzing,  as  a  typical  example  of 
direct  conversion  heat  engines,  the  thermoelectric  generator.  The  effects  of  heat-leak  and  finite-rate  heat 
transfer  on  the  performance  of  the  generator,  which  is  modeled  as  a  Carnot-like  engine  with  internal 
irreversibilities,  are  studied.  Even  though  both  methods  lead  to  the  same  conclusions,  the  entropy 
generation  minimization  method,  when  applied  for  such  an  inherently  irreversible  device,  is  shown  to  be 
less  straight  forward  than  the  power  maximization  technique  requiring  careful  accounting  of  the 
different  sources  of  irreversibility.  Moreover,  the  entropy  generation  rate  vs.  efficiency  behavior  of  the 
generator  reveals  that  the  efficiency  at  minimum  entropy  generation  rate  and  the  maximum  efficiency  are 
distinct.  ©  2000  Elsevier  Science  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  basic  tools  for  studying  heat  engine  performance  are  based  on  aspects  of  the  first  and 
second  laws  of  thermodynamics.  In  fact,  in  terms  of  the  simplest  possible  considerations, 
complicated  thermal  systems  can  easily  be  assessed.  While  detailed  predictions  of  the  problem 
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Nomenclature 

/  Electric  current  (A) 

K,  Kc,  A'h  Heat  transfer  coefficients  (W/K) 
m  Load  ratio 

Q\  Heat  entering  engine  (W) 

Q2  Heat  leaving  engine  (W) 

Qc  Heat  entering  heat  sink  (W) 

Qh  Heat  leaving  heat  source  (W) 

Qe  Excess  heat  bypassing  generator  (W) 

R  Electrical  resistance  (Q) 

Sgen  Entropy  generation  rate  (W/K) 

T\  Temperature  of  hot  junction  (K) 

73  Temperature  of  cold  junction  (K) 

7c  Temperature  of  heat  sink  (K) 

Th  Temperature  of  heat  source  (K) 

W  Power  output  (W) 

Z  Figure  of  merit  (K_1) 

a  Seebeck  coefficient  (V/K) 

r/  Efficiency  of  heat  engine 

rjc  Carnot  efficiency 

?ymax  Maximum  efficiency 

r\mp  Efficiency  at  maximum  power 


considered  are  not  possible,  a  rather  simplified  view  of  the  system  fulfilment  can  be  obtained 
which  may  be  of  value  to  a  designer  trying  to  choose  among  several  suggested  arrangements 
before  going  into  detailed  analysis. 

The  performance  of  heat  engines  has  been  the  subject  of  many  studies  using  what  appear  to 
be  different  approaches,  while  in  fact,  all  these  methods  are  similar  and  can  be  classified  to  fall 
under  the  Power  Maximization  (PM)  technique.  Bejan  [1]  recently  reviewed  the  available 
procedures  and  suggested  a  new  approach  called  “entropy  generation  minimization”  (EGM). 
The  method,  a  combination  of  thermodynamics,  heat  transfer  and  fluid  flow  essentially  unifies 
what  is  known  in  engineering  as  “thermodynamic  optimization”  and  in  physics  as  “finite-time 
thermodynamics”.  Moreover,  in  the  widely  used  PM  method,  the  heat  engine  performance  is 
optimized  by  maximizing  power  generation,  while  in  the  EGM  method,  the  same  task  is 
accomplished  by  minimizing  entropy  generation. 

Several  workers  used  the  concept  of  endoreversibility  on  a  heat  engine  to  predict 
performance.  Among  others,  Curzon  and  Ahlborn  [2]  and  DeVos  [3]  applied  finite-time 
thermodynamics  to  analyze  heat  engines  operating  at  maximum  power.  Gordon  and  Zarmi  [4] 
and  DeVos  and  Flater  [5]  modeled  the  earth  and  its  envelope  as  a  Carnot-like  engine  with  heat 
input  being  solar  radiation  and  work  output  the  wind  energy  generated.  Nuwayhid  and 
Moukalled  [6]  looked  further  into  the  addition  of  a  heat-leak  term  and  studied  the  effect  of  a 
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planet  thermal  conductance  on  the  conversion  efficiency  of  solar  energy  into  wind  energy. 
Moukalled  et  al.  [7]  developed  a  generalized  model  in  which  a  heat-leak  term  is  added  to 
account  for  friction  as  well  as  external  thermal  losses  from  the  engine  to  the  surroundings  in 
addition  to  heat  exchange  between  the  engine,  heat  reservoirs  and  the  surroundings  via 
combined  modes  of  heat  transfer. 

The  thermoelectric  generator  represents  a  typical  example  of  an  inherently  irreversible,  direct 
conversion  heat  engine  and  entails  interesting  electrical  as  well  as  thermal  phenomena. 
Optimization  studies  have  usually  concentrated  on  power  assessment  and  rarely  on  entropy 
tracking  criteria.  An  early  paper  on  the  second  law  evaluation  of  the  thermoelectric  generator 
was  that  of  Lampinen  [8].  Lampinen’s  treatment  was  of  value  since  it  included  the  Thomson 
heat,  but  it  did  not  address  finite-rate  heat  transfer,  nor  did  it  account  for  the  power 
production  dependence  on  the  load.  Gordon  [9],  using  the  thermoelectric  generator  as  an 
example,  clearly  demonstrated  that  real  heat  engines  with  sources  of  irreversibility  that  include 
friction  and  heat-leak  exhibit  fundamentally  different  power  vs.  efficiency  curves  than  those 
predicted  neglecting  friction  and  heat-leak.  Nuwayhid  and  Moukalled  [10]  extended  Gordon’s 
view  of  the  thermoelectric  generator  and  recently  reported  on  the  performance  assuming 
different  thermodynamic  models  and  conditions.  The  present  work  attempts  to  look  into  the 
entropy  generation  rate  suffered  by  thermoelectric  power  devices  and  to  compare  the  PM 
method  vs.  the  EGM  method.  The  ultimate  goal  of  this  work  is,  however,  to  further  promote 
utilization  of  the  thermoelectric  effect. 

In  the  remainder  of  this  paper,  the  first  and  second  laws  of  thermodynamics  are  first 
reviewed.  Then,  the  traditional  treatment  of  a  thermoelectric  generator,  as  an  internally 
irreversible  heat  engine,  using  the  power  maximization  technique  is  elaborated.  This  is  followed 
by  a  detailed  study  of  the  generator  using  the  EGM  method.  The  EGM  analysis  is  performed 
on  three  different  models  of  the  thermoelectric  generator.  In  the  first,  the  generator  is  taken  to 
be  internally  irreversible  and  externally  reversible  while  in  the  second,  the  opposite  is  assumed. 
The  third  model  is  the  most  general  one  in  which  the  generator  is  considered  to  be  internally 
as  well  as  externally  irreversible.  For  all  three  models,  a  comparison  between  the  EGM  and 
PM  analyses  is  performed. 


2.  The  laws  of  thermodynamics 

In  general,  the  first  and  second  laws  for  any  selected  “system”  may  be  written  as  follows: 

^7  =  +  £z)  -  Yjk{h  +  Kr  +  gz)  +  Q  —  W  (1) 

01  In  V  Z  7  Out  v  z  7 

and 

■  dS  ^  ^  Q 

Sgen  =  -  Yjns  +  Yjns  -  —  (2) 

dt  o^t  ro 


where  the  terms  are  those  found  in  any  standard  thermodynamics  textbook.  At  steady  state 
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(time  derivatives  are  zero)  and  for  a  closed  system  (in/out  flow  terms  zero),  the  laws  become 
(assuming  there  may  be  more  than  one  thermal  exchange  between  the  system  and  surroundings 
—  “N”  such  contacts): 

N 

W=  Y^Qi  (3) 

i=\ 

and 

The  above  two  equations  are  applicable  to  any  “system”  residing  within  a  temperature 
difference.  Bejan  [1]  very  appropriately  referred  to  a  heat  engine  operating  in  a  finite 
temperature  difference  field  as  simply  being  an  “insulation”.  This  wonderfully  simple  concept 
allows  very  useful  observations  to  be  made.  In  fact,  Bejan  [1],  through  the  work  of  Chambadal 
[11]  and  Curzon  and  Ahlborn  [2],  demonstrates  clearly  that  the  maximum  power  and  minimum 
entropy  generation  rate  points  are  equivalent.  The  study  of  the  entropy  buildup  suffered  by  a 
heat  engine  provides  a  good  demonstration  of  the  useful  power  one  can  expect  from  such 
devices.  This  is  exploited  in  this  work  in  an  attempt  to  optimize  the  performance  of  that 
delightfully  simple  under-used  power  generating  device:  the  thermoelectric  generator. 


3.  The  PM  method  applied  to  the  thermoelectric  generator 

The  thermoelectric  effect  was  first  reported  by  Seebeck  in  1821  when  he  observed  that  upon 
the  joining  of  two  dissimilar  metal  wires  and  the  heating  of  the  junction,  an  electromotive  force 
developed.  Peltier  in  1834  found  that  the  crossing  of  an  electric  current  across  a  junction  of 
two  dissimilar  conductors  produced  a  heating  or  cooling  effect  depending  on  the  current 
direction.  In  1856,  Thomson  (Lord  Kelvin)  showed  that  there  exists  a  unique  relation  between 
the  Seebeck  and  Peltier  effects  and  while  so  doing  he  actually  discovered  a  totally  new  effect  — 
the  Thomson  effect.  In  this  work  the  Seebeck  and  Peltier  effects  will  be  utilized  in  the 
formulation,  while  the  Thomson  effect  will  be  neglected  based  on  the  assumption  that  the 
properties  are  constant  at  the  average  temperature. 

The  thermoelectric  generator,  having  a  low  thermal  to  electric  conversion  efficiency  has  not 
managed  to  make  its  place  in  the  commercial  market.  Only  in  temperature  measurement  has 
there  been  a  significant  application  of  the  thermoelectric  phenomenon.  Advances  in  the  field 
have  been  very  slow  and  limited  by  material  considerations.  However,  with  the  introduction  of 
semiconductors  and  the  resultant  relative  rise  in  efficiency,  renewed  interest  was  stirred  but 
today  the  practical  applications  still  remain  few.  Current  applications  are  in  remote  low  level 
power  generation  using  long-lived  heat  sources  such  as  radioisotopes,  and  in  small  portable 
refrigeration  systems  running  on  low  power  inputs  such  as  car  batteries.  As  a  result  of  all  this, 
the  few  existing  manufacturers  of  thermoelectric  devices  are  today  finding  it  hard  to  survive 
and  make  a  profit. 
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The  basics  of  thermoelectric  power  generation  can  be  found  in  several  available  texts,  for 
example,  Soo  [12].  It  is  of  interest  to  note  that  in  the  thermoelectric  generator,  the  heat-leak 
irreversibility  is  a  necessary  part  of  the  engine,  since  it  is  a  manifestation  of  the  thermal 
conductance  of  the  generator  legs  themselves  and,  therefore,  appears  in  the  modeling  directly. 
Moreover,  the  thermal  conductance  is  inversely  proportional  to  the  electrical  resistance  and, 
therefore,  the  two  have  conflicting  effects  —  the  system  is,  thus,  one  that  could  be  described  as 
“inherently  irreversible”.  A  simplified  Schematic  of  a  thermoelectric  generator  is  depicted  in 
Fig.  1. 

In  most  scoping  studies  of  thermoelectric  generation,  the  usual  practice  has  been  to  neglect 
external  irreversibilities  and  to  consider  only  internal  effects  that  are  related  to  the  selected 
material.  However,  a  complete  thermodynamic  study  may  consider,  in  addition  to  the 
internally  irreversible  case,  both  externally  irreversible-internally  reversible  and  fully 
irreversible  cases.  Since  the  intention  in  this  section  is  to  give  an  overview  of  the  analysis  using 
the  PM  technique,  the  fully  irreversible  case  is  not  considered  and  attention  is  focused  on  the 
standard  internally  irreversible -externally  reversible  and  internally  reversible-externally 
irreversible  cases. 

With  external  irreversibilities  neglected,  the  heat  source  and  heat  sink  temperatures  become 
identical  to  the  hot  junction  and  cold  junction  temperatures,  respectively.  The  resulting 
equations  for  power  and  efficiency  are  usually  written  as: 


Heat  Source,  TH 


Fig.  1 .  A  schematic  of  a  basic  thermoelectric  heat  engine. 


896 


R.Y.  Nuwayhid  et  al.  /  Energy  Conversion  &  Management  41  (2000)  891-914 


W=  a I(Th  -  7c)  -  I2R, 


(5) 


n  =  nc — 

i  + 


AT  i  (1+m)2 
27fi  +  ZTh 


(6) 


Where  m  =  Ry/R,  Z  =  y2/ RK,  R,  K,  a  and  rjc  are  the  load  ratio,  the  figure  of  merit,  the  total 
electrical  resistance,  the  total  thermal  conductance,  the  Seebeck  coefficient  and  the  Carnot 
efficiency,  respectively. 

Eq.  (6)  is  the  usual  one  found  in  textbooks  and  the  power  vs.  efficiency  behavior  is  the 
expected  “loop”  shape  with  both  efficiency  and  power  rising  from  zero  and  returning  to  zero 
after  passing  through  distinct  maximum  power  and  maximum  efficiency  points  below  the 
Carnot  limit.  The  loss  of  power  is  attributed  to  heat-leak  at  open  circuit  (slow  operation;  no 
recovered  output)  and  to  Joule  heat  loss  at  short  circuit  (fast  operation;  friction).  This 
observation  was  first  made  by  Gordon  [9]. 

If  the  generator  is  made  more  and  more  reversible  (say  by  making  R  smaller  and  smaller 
while  holding  K  artificially  constant  for  simplicity),  the  graph  starts  to  approach  the  “dome” 
shape  where  the  maximum  efficiency  is  the  Carnot  efficiency  and  the  maximum  power  and 
maximum  efficiency  points  become  increasingly  separated  (Fig.  2  shows  this  behavior).  The 
useful  power  can  be  expressed  in  terms  of  m  as  follows: 


Fig.  2.  Typical  effect  of  increasing  internal  irreversibility  (in  the  direction  of  the  arrow)  on  the  performance  of 
thermoelectric  generator  (7h  =  600  K,  t]c  —  0.5). 
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W=  ZK(Th  -  Tcf  m 

(1  +mf 

which  clearly  exhibits  a  maximum  at  m  —  1,  as  matched  load  conditions  require. 
From  Eq.  (6),  with  m  =  1,  the  efficiency  at  maximum  power  is  given  by: 

2 

Imp  =  nc - g 
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(7) 


(8) 


The  previous  figure  (Fig.  2)  shows  that  as  the  device  becomes  more  reversible,  the  maximum 
power  line  shifts  to  the  right  and  approaches  a  limiting  value  which  is  only  dependent  on  the 
Carnot  efficiency.  The  limit  is  obtainable  by  requiring  Z— >oo,  giving: 

=  (9) 

The  above  equation  raises  an  interesting  question  about  the  “maximum”  efficiency  at 
maximum  powers  which  appears  to  be  limited  to  2/3! 

The  maximum  efficiency  is,  however,  distinct  from  the  efficiency  at  maximum  power  and  is 
best  written  as: 


Load  Ratio  (m) 

Fig.  3.  Renewable  efficiency  compared  to  exhaustable  efficiency  for  an  internally  irreversible  generator  (7h  =  600  K, 
r\c  =  0.5,  Z  =  0.001  and  K  —  0.1).  For  the  renewable  case  (solid  lines),  the  inexhaustible  heat  supply  Q  decreases 
from  bottom  to  top  as:  90,  50  and  39. 
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(10) 


where  f  is  the  average  temperature  along  the  generator.  From  Eq.  (10),  it  can  easily  be 
inferred  that  the  load  ratio  giving  maximum  efficiency  is  not  unity  but  v  1  +  ZT  a  value  greater 
than  unity.  Of  course,  as  Z— >oo,  the  Carnot  limit  is  approached.  Further,  it  is  obvious  that  the 
maximum  efficiency  is  distinct  (and  larger)  than  the  efficiency  at  maximum  power  (see  Fig.  2). 

It  may  be  of  interest  to  compare  the  exhaustible  efficiency  given  by  Eq.  (6)  with  the  so-called 
“renewable”  efficiency  which  may  be  given  by: 


mKZT&il 

(1  +  m)2Q 


(11) 


where  Q  is  an  “inexhaustible”  or  “renewable”  heat  source.  While  a  derivative  of  this  with 
respect  to  m  will  show  this  to  maximize  at  m  —  1 ,  the  exhaustible  efficiency  maximizes  when 
m  =  V 1  +  ZT.  Fig.  3  shows  the  exhaustible  efficiency  contrasted  with  renewable  efficiencies  at 
three  different  values  of  Q.  The  maximum  renewable  efficiency  which  occurs  when  Q  = 
K(Th  -  7c)(l  +  ZTC)  is  given  by: 


>r0fn  -Tc) 

1,1  (l+m)2(l+zrc)’ 


(12) 


and  is  shown  to  be  less  than  the  exhaustible  efficiency. 

With  no  external  irreversibilities  considered,  improvement  in  performance  can  proceed 
through  material  selection.  In  fact,  there  has  been  some  work  on  thermoelectric  material 
selection,  through  either  higher  Seebeck  coefficients  or  by  optimizing  the  geometrical 
configuration  of  the  device,  leading  to  larger  Z  and,  thus,  greater  efficiencies. 

The  effect  of  external  irreversibilities,  on  the  other  hand,  is  the  same  for  any  internally 
reversible  (i.e.  endoreversible)  heat  engine  and  has  been  studied  intensely  in  the  literature  [1-6]. 
Although  this  model  is  highly  unrealistic  for  the  present  application,  it  nevertheless  serves  as  a 
link  among  models.  For  the  thermoelectric  generator,  both  R  and  K  are  required  to 
simultaneously  vanish  (in  reality  they  are  inversely  proportional),  and  the  resulting  energy 
balances  become  very  simple.  The  power  vs.  efficiency  relation  is  found  to  be: 

w-  (13) 

(l-rtfsTH  +  O-rf'-'W] 

where  i  is  the  heat  input/output  power  mode  (1  for  conduction/convection  or  4  for  radiation). 
The  normalized  power  (W*  =  W/KT‘H  )  vs.  efficiency  curves,  with  the  Carnot  efficiency  and  the 
heat  transfer  mode  as  parameters,  have  the  expected  “domed”  shape  with  the  maximum 
efficiency  terminating  at  the  Carnot  limit  (by  differentiating  the  above  with  respect  to  K  and 
setting  to  zero  one  indeed  obtains  the  Carnot  limit).  The  efficiency  at  maximum  power  in  this 
“endoreversible”  case,  obtainable  this  time  by  differentiating  the  above  with  respect  to  rj  and 
setting  to  zero,  turns  out  to  be  the  so-called  CNCA  [3]  efficiency  given  by: 
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The  limit  of  this  as  7fi— »oo  is  unity,  compared  to  the  upper  limit  of  2/3  for  the  internally 
irreversible  case.  Moreover,  it  turns  out  that  nearly  below  r\c  =  0.75,  the  endoreversible  (with 
reversibilities  occurring  externally)  and  internally  irreversible  (with  no  external  irreversibilities) 
cases  produce  very  similar  efficiencies  at  maximum  power  as  long  as  Z  is  greater  than  about 
0.1. 

As  a  general  guideline  consider  a  thermoelectric  generator  with  the  following  parameters: 
7fi  =  600  K,  Tc  =  300  K,  a  =  10“4  V/K,  R  =  10“5  Q  and  K  =  0.1  W/K.  In  the  fully  reversible 
limit,  the  Carnot  limit,  the  efficiency  is  obviously  50%.  The  externally  reversible  assumption 
gives  20%  to  24%  efficiency  as  calculated  by  Eqs.  (8)  or  (10)  respectively.  The  endoreversible 
case  (CNCA  efficiency)  gives  30%  for  the  efficiency  at  maximum  power.  All  these  figures  are 
very  optimistic  and  real  generators  seldom  exceed  10%  efficiency. 


4.  The  EGM  method  applied  to  the  thermoelectric  generator 

An  overview  of  the  performance  of  the  thermoelectric  generator  using  the  PM  technique  was 
given  in  the  previous  section.  In  this  section,  the  generator  will  be  investigated  using  the  EGM 
method.  Entropy  generation  in  a  thermoelectric  device  occurs  with  heat  transfer  from  the  heat 
source,  heat  transfer  to  the  surroundings,  heat-leak  from  the  hot  side  to  the  cold  side  and 
dissipative  Joule  heat.  The  irreversibilities  may  occur  internally  or  externally.  The  analysis  will 
proceed  from  simplified  models  to  a  general  one.  In  specific,  three  models  will  be  studied:  (1) 
the  internally  irreversible  and  externally  reversible  model,  (2)  the  internally  reversible  and 
externally  irreversible  model  and  (3)  the  internally  and  externally  irreversible  model.  The  prime 
motive  behind  this  work  is  the  apparent  lack  of  studies  on  entropy  generation  in  thermoelectric 
devices. 


4.1.  Entropy  generation  in  the  internally  irreversible-externally  reversible  model 

The  generator  is  modeled  as  externally  reversible  but  internally  irreversible.  The  heat  transfer 
to  the  generator  (gH)  and  from  the  generator  (Qc)  are: 

Qn  =  a ITh  +  K(Th  -  Tc)  -  I2R/ 2  (15) 

and 

Qc  =  *ITc  +  K(T}i-Tc)  +  I2R/2  (16) 

where  the  Joule  heat  has  been  split  evenly  between  the  two  sides  of  the  generator.  This  usually 
implies  that  the  Joule  heating  effect  is  larger  in  magnitude  than  the  heat  conduction  effect 
along  the  generator. 

The  entropy  generation  rate  is,  as  usual,  Sgen  =  — X^=i(f)p  where  the  i  refers  to  the  exchange 
location  (node).  So  that,  for  this  apparently  two  node  case: 
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Sgen  — 


Qh  ,  Qc 
Th  Tc 


(17) 


Applying  the  above  equations  gives: 


1  «  Jii  -j-  Tc 

K(TH-TC)  +  -I2R  c 

L  2  yH  -  2c j 


(18) 


This  result  is  in  conformity  with  Bejan  [1]  in  reference  to  entropy  generation  in  a  general 
temperature  field  with  no  thermoelectric  effect  and,  thus,  no  electric  current  flow.  The 
irreversibility  in  the  thermoelectric  generator  is  clearly  due  to  the  conduction  of  heat  and  to  the 
dissipation  of  Joule  heat  (the  Seebeck/Peltier  effects  do  not  give  rise  to  entropy  creation  since 
they  are  reversible  effects).  Moreover,  the  conduction  heat  given  by  QK  =  K(TH  —  Tc)  generates 
entropy  in  the  amount  gK(^r~^c),  while  the  Joule  heat  (Q,  =  I2R )  generates  entropy  in  the 
amount:Qj(iJg;)- 

The  equation  for  the  total  entropy  generation  rate  can  best  be  written  in  dimensionless  form 
as: 


^  =  (ti/2  +  t-i/2)2  i  ,  j 

K  |_  2(1  +  m)2 


(19) 


where  t  is  the  temperature  ratio  (7~c/7h)  and  Z  and  m  are  as  defined  earlier.  It  is  seen  that  as 
Z  gets  larger,  more  entropy  is  generated  due  to  the  added  effect  of  the  Joulean  heating 
resulting  from  the  induced  current.  Additionally,  as  m  gets  larger,  R  gets  smaller  and  the 
generator  is  less  irreversible  thereby  generating  less  entropy.  There  is  a  problem,  however,  with 
what  has  been  done  so  far!  Strictly  speaking,  what  is  being  assessed  is  the  case  where  the 
irreversibilities  are  only  internal,  and  so,  any  external  sources  of  entropy  generation  are  not 
considered.  However,  to  make  the  study  more  realistic  and  consistent  with  the  traditional 
treatment  of  thermoelectric  performance,  a  little  more  thought  is  required. 

Looking  back  at  the  previous  equation  for  entropy  generation  rate,  it  is  immediately 
apparent  that  it  does  not  exhibit  the  appropriate  minimum  entropy  at  m  =  1 .  This  is  of 
concern,  since  the  maximum  power  and  minimum  entropy  points  ought  to  be  one  and  the 
same  [1].  The  problem  can  be  traced  back  to  neglect  of  the  entropy  generation  due  to  external 
heat  dumping  (Bejan’s  room  to  move  criteria).  Including  this  effect,  one  considers  that  a  source 
of  heat  Q  is  available  at  Tfi  and  a  part  of  this,  gH  enters  the  generator  at  Tfi,  while  another 
part  Qt  is  dumped  externally  towards  the  sink  temperature  Tc.  As  such,  one  may  write  the 
entropy  generation  rate  as  the  sum  of  those  of:  (1)  internal  heat-leak,  (2)  internal  Joule  heating, 
and  (3)  external  heat  dumping.  The  following  is,  thus,  obtained: 


Sgen 


=  K 


(Th  ~  Tc) 
ThTc 


-i2r 


Th  +  Tq 
ThTc 


+  a 


(Th  ~  Tc) 

'■  ThTc 


(20) 


The  above  equation  can  alternatively  be  written  as: 
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V  t  +  t+K  Th+  Tc)  (21) 

After  writing  Qe  =  Q  —  Qn,  replacing  the  current  using  tx.1  —  KZ(TH  —  Tc)/(1  +  m)  and 
performing  some  algebraic  manipulations,  either  of  the  previous  two  equations  can  be 
presented  in  the  following  form: 


a  _  a(Th  -  Tc)  _  mKZ(TH  -  Tc  f 
gen  Q  THTC  (1  +m)2Tc 


(22) 


The  above  equation  can  be  differentiated  with  respect  to  m  to  show  clearly  that  a  minimum 
occurs  at  m  =  1.  Thus,  maximum  power  and  minimum  entropy  are  one  and  the  same  point. 
For  further  clarity,  the  equation  is  shown  plotted  in  Fig.  4  in  the  following  dimensionless  form: 


■*  _  mZTH  K(Th  -  Tc) 

gen  +m)2  'M:  * 


(23) 


where  the  last  term  may  be  attributed  to  the  fraction  of  the  heat  source  that  goes  into  the  heat- 
leak.  The  figure  shows  clearly  that  as  the  heat-leak  increases,  the  entropy  generation  rate 
increases  and  the  useful  power  diminishes.  Additionally,  Fig.  5  is  presented  to  show  how  the 
increase  in  the  thermoelectric  figure  of  merit  improves  performance  (i.e.  reduces  entropy 
generation  rate). 


Load  Ratio  (m) 

Fig.  4.  Dimensionless  entropy  generation  rate  vs.  load  ratio  for  several  values  of  heat  leak  in  the  internally 
irreversible  model  (TK  =  600  K,  Tc  —  300  K  and  Z  —  0.001).  The  heat  leak  increment  (K(Ta  —  Tc)/Q)  increases 
from  top  to  bottom  as:  0.01,  0.1,  0.5,  0.9  and  0.99. 
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Load  Ratio  (rri) 

Fig.  5.  Dimensionless  entropy  generation  rate  vs.  load  ratio  in  the  internally  irreversible  model  for  three  values  of 
thermoelectric  figure-of-merit  Z  increasing  from  top  to  bottom  as:  0.001,  0.002  and  0.005  (Th  =  600  K,  7c  =  300  K 
and  K(Th  -  TdQ  =  0.5). 


Fig.  6.  Equivalent  model  of  the  internally  irreversible/externally  reversible  thermoelectric  generator. 
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Combining  Eqs.  (7)  and  (22),  the  following  linear  relation  is  found: 

.  (T„  -  Tr) 

W  —  Qc  H  CJ-TcSgen  (24) 

7h 

The  equation  shows  clearly  that  of  the  maximum  available  power  (first  term  on  the  right),  an 
amount  equal  to  the  second  term  on  the  right  is  lost. 

Eq.  (24)  can  be  recast  as: 

^  =  (25) 

which  states  that  in  this  model,  the  entropy  generation  rate,  as  schematically  shown  in  Fig.  6, 
is  equivalent  to  an  inexhaustible  heat  source  ( Q )  proceeding  undiminished  through  a 
temperature  difference  (Tfi  —  Tc),  in  parallel  with  a  Carnot  engine  producing  a  power 
(2h  —  2c)-  The  result  is  consistent  with  the  Guoy  Stodola  theorem  on  destruction  of 
availability.  Power  is  entropy  free  and  the  diagram  shows  the  produced  power  dissipated  into 
Tc- 

Next,  a  consideration  of  the  case  of  the  internally  reversible/externally  irreversible 
thermoelectric  generator  may  afford  further  insight. 


4.2.  Entropy  generation  in  the  internally  reversible-externally  irreversible  model 


Here,  irreversibility  is  assumed  to  occur  externally  and  one  has  the  so-called  endoreversible 
case.  This  application  makes  little  real  sense,  since  thermoelectric  devices  are  internally 
irreversible  by  nature.  It  does,  however,  provide  a  link  with  both  the  fully  irreversible  case  and 
the  general  thermal  power  generation  plant  models  [1]. 

The  heat  flow  to  the  hot  junction  from  the  heat  source  and  from  the  cold  junction  to  the 
heat  sink  are: 


qh  =  Kh(Th-T1)  and  QC  =  KC{T2-TC )  (26) 

Additionally,  there  is  assumed  to  be  no  internal  irreversibilities  (R  =  0,  K  =  0)  and  the  heat 
flows  into  the  generator  and  out  of  it  are  conserved  and  given  by: 

2h  =  olIT\  and  Qc  —  aIT2  (27) 

The  above  equations  can  be  solved  simultaneously  to  yield: 


KhTh 
1  a/+  Kh 


(28) 


KCTC 
Kq  —  a/ 


(29) 


By  eliminating  a  I  between  the  two  equations,  the  following  relation  between  7)  and  T2  is 
obtained: 
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KCTC 


Kh  +  Kc  —  Kyi 


m 


(30) 


The  total  entropy  generation  rate  is  equal  to  the  sum  of  the  rates  of  entropy  generated  by 
finite-rate  heat  transfer  at  the  junctions  and  that  generated  by  the  excess  heat  dumping  (Qc) 
due  to  the  inexhaustible  heat  source  Qe(Q—Qe  +  QH )  and  is  given  by: 


Sgen  —  Q  H  | 


'Th-TA 


This  can  also  be  written  as: 

((Tc  -  T0(TH  -  TO 


S’  en  =  Kf 


TxTc 


(31) 


(32) 


(33) 


Replacing  T2  by  its  equivalent  expression  given  by  Eq.  (30),  the  entropy  generation  rate 
becomes: 


JTU-TC\  Kh(Th-T0{Kc(Tc-T0  +  Kh(Th-T0) 
V  ThTc  )  TOKhT  +  KCT\  -  KhTh) 


(34) 


Taking  dSgen/dT\  and  setting  to  zero,  the  optimum  hot  side  temperature  that  minimizes  the 
entropy  generation  rate  is  found  to  be: 


T i  _  opt  — 


Kq  +  Kh 


(35) 


Substitution  into  Eqs.  (30)  and  (33)  results,  respectively,  in: 

T ,  _  Kh^/ThTq  +  KCTC 

r2’°pt“  K^+Kh 


(36) 


and 


*Sgen,  min  —  Q  | 


Th-Tc\  Kh  (V7^-V7d)2 
ThTc  )  +  Kh  +  Kc  Tc 


(37) 


Based  on  the  obtained  values  for  T\  and  T2  that  minimize  the  rate  of  entropy  creation,  the 
output  of  the  heat  engine  is  given  by: 


(38) 
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and  the  efficiency  is 

_  i  _  Kqy/  Th  Tc  +  Kc  Tq 
'?opt  “  kcVt^  +  khth 

It  can  be  seen  that  as: 

KC—>00,  Tyopt^-^/TcTu,  T2.  opt^Tc  and  ^gen,  min— 7^0 Th  C  ) 
On  the  other  hand,  as: 


(39) 


(Tn-Tc\  (V7k-V7b)2 


Similar  conclusions  could  be  made  regarding  the  optimum  efficiency. 

Comparing  the  optimum  temperatures  and  power  obtained  here  via  the  EGM  method  with 
those  obtained  in  DeVos  [3]  using  the  PM  technique,  they  are  found  to  be  equal,  confirming 
that  minimum  entropy  generation  rate  and  maximum  power  output  conditions  are  coincident. 


Fig.  7.  Gc 
generator. 
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4.3.  Entropy  generation  in  the  internally  and  externally  irreversible  model 

In  the  previous  two  Carnot-like  engines  used  to  model  the  thermoelectric  generator,  either 
the  internal  or  the  external  irreversibility  was  neglected.  In  this  section,  a  general  model  that 
accounts  for  both  irreversibilities  is  sought.  For  that  purpose,  the  Carnot-like  engine’s  model 
developed  by  Moukalled  et  al.  [13]  is  modified  to  account  for  internal  irreversibilities  within  the 
engine  and  used  here  to  analyze,  following  the  EGM  approach,  the  thermoelectric  generator.  A 
Schematic  of  the  modified  model  is  depicted  in  Fig.  7.  As  can  be  seen,  in  addition  to  the 
internal  irreversibilities  and  heat-leak,  the  model  allows  for  finite-rate  heat  input  and  output 
with  heat  transfer  obeying  Newton’s  law. 

Referring  to  Fig.  7,  the  heat  balance  equations  at  the  junctions  are  given  by: 

Ah(7h  -Tx)  =  ocIT1-I^  +  K(T}  -  T2)  (40) 


i2r 


(41) 


where  the  terms  to  the  left  of  the  equality  are  QH  and  Qc,  while  those  to  the  right  are  called 
Qx  and  Q2  and  are  the  internal  heat  flows  (note,  of  course,  that  the  first  law  requires  QH  =  Q\ 
and  Qc  =  Q2). 

Solving  Eqs.  (40)  and  (41)  simultaneously,  the  following  equations  for  7j  and  T2,  as  a 
function  of  the  current  /,  are  obtained: 


T\  = 


2KKCTC  +  2 KKkTk  +  2KcKhTh  -  2KHTHI+  (2 KR  +  KCR)I2  -  Roe/3 
2(KKC  +  KKh  +  KcKh  +  IKcy.  -  IKHy.  -  /2a2) 


(42) 


2KKC  Tq  +  2KKhTh  +  2KcKhTc  +  2KCTCI+  (2  KR  +  KHR)I2  +  Ra/3 
2 (KKC  +  KKh  +  KcKh  +  lKcy.  -  IKHy.  -  Py2) 


(43) 


The  entropy  generation  rate  may  now  be  found  by  noting  that  this  case  is  a  multi-node  case. 
Thus,  in  addition  to  the  finite-rate  temperature  difference  irreversibilities  (first  two  terms 
below),  there  are  the  internal  irreversibilities  due  to  the  heat-leak  and  Joule-heat  (the  sum  of 
the  third  and  fourth  terms)  in  addition  to  the  irreversibility  due  to  external  heat  dumping 
across  the  imposed  temperature  field.  Thus: 


Sgen  -  +  GC(%^)  -  ^  +  6c| 


7Wc\ 

TyyTc  ) 


(44) 


As  stated  before,  the  first  law  requires  QH  and  Qc  to  be  equal  to  Q\  and  Q2  respectively.  In 
addition,  as  in  the  previous  cases,  the  external  dumping  is  written  as  Qe  —  QH  —  Qc  where  QH 
is  a  large  heat  supply  available  at  TH.  After  some  manipulation,  the  resulting  entropy 
generation  rate  becomes: 
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where  the  last  quotient  of  the  second  term  can  be  viewed  as  the  fraction  of  the  heat  supply 
that  goes  into  the  heat-leak.  If  a  derivative  with  respect  to  m  is  now  taken,  one  obtains  mopt  = 
1  as  expected. 

If,  on  the  other  hand,  one  writes: 


Sgen  —  Q 


( 


Tn-Tc 

ThTc 


Kh(Th  -  Tx)  -  Kc(T2  -  7c) 
Tc 


(49) 


and  if  the  values  for  T\  and  A?,  as  given  by  Eqs.  (42)  and  (43),  are  substituted  into  Eq.  (48), 
the  following  equation  for  the  entropy  generation  rate  is  obtained: 


(Kc  -  Ka)Rcd3  +  2ACAH(AC  -  T„)y. 
2AC(A(AC  +  Ah)  +  (Ac  -  a/)(AH  +  a/)) 


2(KcTca2  +  AHrHo i2  +  KcKnR  +  K(KC  +  KH)R)l2 
2Tc{K(Kc  +  Ah)  +  (Ac  -  «7)(AH  +  a I)) 


(50) 


Now,  if  the  derivative  with  respect  to  current  /  of  the  entropy  generation  rate  is  set  to  zero,  the 
current  minimizing  the  entropy  generation  rate  emerges.  The  resultant  relation  for  this  current 
is  very  long  and  it  would  be  of  little  value  to  report  it  here.  Rather,  if  the  reasonable 
assumption  Ac  =  AH  is  used,  then  the  analytical  expression  for  the  optimum  current  becomes 
manageable  to  present  and  is  found  to  be: 


^  yj (2A  +  AH) { (2A  +  AH)[(TC  +  TH)a2  +  2KR  +  KhR]2-Kh(Tc  -  7h)V  } 


(Tc  -  Ah) a3 

(2A+  AH)  [(Ac  +  AH)a2  +  2KR  +  AHi?] 
(Ac  -  AH) a3 


(51) 


Using  /opt,  the  following  expressions  for  Ty  and  entropy  generation  rate  are  obtained: 

T\,  opt  = 

{(Ac  -  AH)a2  +  4KR  +2KHR}^(2K+KH){(2K+Kli)[(Tc+Tii)<x2  +2KR  +  AhR]2-Ah(Ac  -  AH)2a4} 
4(2A+Ah)(Ac-  Ah)cc4 


16 A2 A 3  +  4AA2{(3Ac  +  AH)a2  +  6AhA}  +  2A{(A2[  -  A2)a4  +  2AHA(3AC  +  AH)a3  +  6A2  A2} 
4(2A+Ah)(Ac- AH)a4 


Ah{2Ah(Ah  -  Ac)a4  +  AHA(3AC  +  AH)a2  +  2A2  A2} 
4(2A+Ah)(Ac- AH)a4 


(52) 
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7 Tu  -  Tc\  KU(2K+  Kh){(Tc  +  7h)oc2  +  2KR  +  KHR ) 

V  ThTc  )  2(2K+  KU)TC^ 


+ 


KHy/(2K+  Kh)(2K+  Kh)(Tc  +  7h)oc2  +  2KR  +  KHRf-KH(Tc  -  TH)2a4) 
2(2  K  +  KH)Tcy.2 


(53) 


The  equation  for  the  optimum  T2 ,  opt  can  also  be  found  but  will  not  be  reported  here  for 
brevity.  Moreover,  the  value  for  T\,  as  given  by  Eq.  (51),  that  minimizes  entropy  is  exactly  the 
value  that  maximizes  power  which  is  derived  in  Ref.  [10]  using  the  PM  technique.  Therefore, 
once  again  and  for  the  general  model,  the  maximum  power  output  is  associated  with  the 
minimum  entropy  generation  rate. 

Using  the  above  derived  equation,  the  effects  of  the  various  parameters  on  entropy 
generation,  power  output  and  efficiency  of  the  generator  are  assessed  and  results  are  displayed 
in  Figs.  9-12.  By  augmenting  the  heat-leak  through  increasing  the  value  of  the  heat  transfer 
coefficient  K  while  holding  all  other  parameters  constant  (Fig.  9),  both  efficiency  and  power 
outpout  decrease,  as  expected,  and  the  rate  of  entropy  generated  increases.  Moroever,  it  is 
obvious  that  the  optimum  current  that  minimizes  the  entropy  generation  rate  maximizes  power 


Fig.  12.  The  effects  of  varying  K  on  the  entropy  generation  rate  vs.  efficiency  in  the  fully  irreversible  model.  Arrow 
shows  K  to  increase  as:  0.001,  0.01  and  0.1  W/K.  ( Tn  =  600  K,  Tc  =  300  K,  a  =  10“4,  R  =  10‘5  and  Kn  =  Kc=  1 .) 
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output.  Furthermore,  the  efficiency  at  minimum  entropy  generation  rate(i.e.  maximum  power 
output)  and  the  maximum  efficiency  of  the  engine  are  distinct. 

The  effects  of  varying  the  Seebeck  coefficient  are  depicted  in  Fig.  10.  As  a  increases,  the 
entropy  generation  rate  decreases  and  the  power  and  efficiency  increase.  This  is  expected  since 
an  increase  in  the  Seebeck  effect  results  in  augmenting  the  “figure  of  merit”  Z  (Z  =  a,2/RK). 
Again,  the  maximum  efficiency  and  minimum  entropy  generation  points  are  not  coincident. 

The  effects  of  finite-rate  heat  transfer  on  the  variations  of  entropy  generation,  power,  and 
efficiency  are  depicted  in  Fig.  11.  As  the  resistance  to  heat  transfer  increases  (i.e.  Kc  decreases), 
the  entropy  generation  rate  increases,  and  the  power  output,  efficiency  at  maximum  power,  and 
maximum  efficiency  of  the  generator  diminish.  This  behaviour  is  anticipated  because  decreasing 
the  thermal  resistance  to  finite-rate  heat  transfer  (by  increasing  Afi  and/or  Kc)  reduces  the 
irreversibilty  of  the  heat  engine  by  forcing  the  hot  and  cold  junction  temperatures  to  become 
closer  to  the  heat  source  and  heat  sink  temperatures  (e.g.  in  this  case  T2  becomes  closer  to  7c). 
As  in  the  previous  plots,  the  minimum  entropy  generation  is  associated  with  the  maximum 
power  output  whereas  the  maximum  efficiency  and  the  efficiency  at  minimum  entropy 
generation  rate  are  distinct. 

Finally,  the  entropy  generation  rate  vs.  efficiency  curves  for  different  heat-leak  values  are 
presented  in  Fig.  12.  As  shown,  the  minimum  entropy  point  and  the  maximum  efficiency  point 
are  distinct  and  by  increasing  the  value  of  K  and  thus  irreversibilty,  they  become  increasingly 
separated.  This  plot  can  be  seen  to  have  the  shape  of  an  inverted  power  vs.  efficiency  curve. 
The  starting  point,  however,  is  not  zero,  as  for  power,  but  is  equal  to  the  entropy  generation 
rate  due  to  the  heat  traverse  of  the  available  temperature  field. 


5.  Concluding  remarks 

A  comparative  assessment  of  the  performance  of  the  thermoelectric  generator  using  the 
EGM  method  and  the  PM  technique  was  presented.  The  effects  of  heat-leak,  internal 
dissipation,  and  finite-rate  heat  transfer  were  investigated.  The  detrimental  effect  of  these 
variables  on  the  performance  of  the  generator  was  demonstrated.  A  tab  on  the  entropy 
generation  rate  suffered  by  the  device  was  kept  and  each  case  was  found  to  have  an  equivalent 
model  consistent  with  the  Guoy-Stodola  theorem.  The  PM  and  EGM  methods  were  shown  to 
lead  to  the  same  conclusions  (i.e.  minimum  entropy  generation  rate  and  maximum  power 
output  conditions  are  coincident)  with  the  EGM  method  being  less  straightforward  than  the 
PM  technique,  requiring  careful  accounting  of  the  different  sources  of  irreversibility.  The  study 
provided  useful  thermodynamic  reasoning  as  applied  to  a  low  maintenance  device  which  is 
under-used  in  a  world  of  bulky  pollution  generating  monsters. 
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